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Abstract 

The Dirichlet Laplacian in a curved three-dimensional tube built along a spatial (bounded or un- 
bounded) curve is investigated in the limit when the uniform cross-section of the tube diminishes. 
Both deformations due to bending and twisting of the tube are considered. We show that the Lapla- 
cian converges in a norm-resolvent sense to the well known one-dimensional Schrodinger operator 

Q, whose potential is expressed in terms of the curvature of the reference curve, the twisting angle and 

a constant measuring the asymmetry of the cross-section. Contrary to previous results, we allow the 
reference curves to have non-continuous and possibly vanishing curvature. For such curves, the dis- 
tinguished Frenet frame standardly used to define the tube need not exist and, moreover, the known 
approaches to prove the result for unbounded tubes do not work. Our main ideas how to establish 
the norm-resolvent convergence under the minimal regularity assumptions are to use an alternative 

i i frame defined by a parallel transport along the curve and a refined smoothing of the curvature via 

the Steklov approximation. 
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1 Introduction 

This paper is concerned with the singular operator limit for the Dirichlet Laplacian in a three-dimensional 
non-self- intersecting curved tube (c/ Figure nj when its two-dimensional cross-section shrinks to a point. 
The tube f2 e is constructed by translating and rotating the cross-section along a spatial curve T and the 
limit is realized by homothctically scaling a fixed cross-section w by a small positive number e. Without 
loss of generality, we assume that the curve is given by its arc-length parameterization r : / — > M 3 , where 
the open interval I C 1 is allowed to be arbitrary: finite, infinite or semi-infinite. Geometrically, O e 
collapses to T as e — > 0. We are interested in how and when the three-dimensional Dirichlet Laplacian 
— A D e can be approximated by a one-dimensional operator H c g on the curve. 




Figure f : The geometry of a quantum waveguide. Twisting and bending are demonstrated on the left 
and right part of the figure, respectively. 



o Since we deal with unbounded operators, the convergence of -A D ' to H c ff is understood through 



o 



We start with some more or less obvious observations. 

Since we deal with unbounded or. 
a convergence of their resolvents. 

The Dirichlet boundary conditions imply that the spectrum of — A B e explodes as e — ¥ 0. It is just 
because the first eigenvalue of the Dirichlet Laplacian in the scaled cross-section ecu :— {at \ t e uj} 
equals e~ 2 Ei, where E\ is the first eigenvalue of the Dirichlet Laplacian in the fixed cross-section u. 
Hence, a normalization -A D ' — e~ 2 E\ is in order to get a non-trivial limit. 

Finally, since the configuration spaces Q e and T have different dimensions, a suitable identification 
of respective Hilbert spaces of —A D e and H e & is required. This is achieved by using a unitary 
transform that identifies L 2 (fl e ) with L 2 (I x lo) and by considering H c g as acting on the subspace 
of L 2 (I x uj) spanned by functions of the form ip ® J\ on I x ui, where J\ denotes the positive 
normalized eigenfunction of — A£, corresponding to E\. 

Taking these remarks into account, we can write the convergence result as follows: 



o 
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A £ _ £ ~2 El f^ H F g ;= _ A T d _ - + Cij {§p T) 2 _ (L1) 



Here K :— |f| and r := n~ 2 det(f ,f , r) denote respectively the curvature and torsion of T, Op is an 
angle function defining the rotation of ew with respect to the Frenet frame of Y and C w := ||9 Q ^i||i^(a))) 
with d a denoting the angular derivative in M 2 . 



The convergence ( 1.1 ) can be employed as a way to approximate the three-dimensional dynamics of an 
electron constrained to a curved quantum waveguide by the effective one-dimensional Hamiltonian H^ s 
on the reference curve. The Dirichlet Laplacian on the interval / represents the kinetic energy of the free 
motion on the reference curve (indeed, — A£ is unitarily equivalent to the Laplace-Beltrami operator 
on r). The additional potential of H^ s clearly consists of two competing terms: the negative one 
induced by curvature and the positive one due to torsion. They respectively represent the opposite 
effects of bending and twisting in quantum waveguides, cf |17j . 



1.1 Known results and why we write this paper 



The result (1.1) is well known, it has been established in various settings and with different methods 



during the last two decades. As the first rigorous result, let us mention the classical paper [S] of Duclos 
and Exner, where the norm-resolvent convergence of (1.1 1 is proved under quite restrictive hypotheses 
r <G C 4 and uj being a disc (so that C u = 0). More precise results about the limit (for instance, uniform 
convergence of eigenfunctions) in arbitrary dimensions are established by Freitas and Krejcifik in |llj . 
but the cross-section is still assumed to be rotated along T in such a way that 8p = t, so there is no 
effect of twisting. 

The presence of the additional potential term due to twisting in H^ s was observed for the first time 
by Bouchitte, Mascarenhas and Trabucho [J]. Contrary to the previous works where operator techniques 
are used, the authors of [1] use an alternative method of Gamma-convergence, which provides just a 
strong-resolvent convergence of (1.1 1 but, on the other hand, enables them to weaken the regularity 
hypothesis to T e C 3 . De Oliveira [gj extended the results of [4] to unbounded tubes and established a 
norm-resolvent convergence in the bounded case (see also 0IE]). 

Finally, let us mention the series of recent papers [20, 23, 24], where the singular limit of the type ( 1.1 ) 
is attacked by the methods of adiabatic perturbation theory. In fact, the general setting of shrinking 
tubular neighbourhoods of (infinitely smooth) submanifolds of Ricmannian manifolds is considered in 
these works and the results can be interpreted as a rigorous quantization procedure on the submanifolds. 



After having provided an extensive literature on the limit ( 1.1 ), a question arises why we still consider 



the problem in the present paper. In fact, the issue we would like to address here is about the optimal 
regularity conditions under which the effective approximation (1.1) holds. We are motivated by the fact 



that the known existing results mentioned above do not cover physically interesting curves with merely 
continuous or even discontinuous curvature ( cf Figure v3 . 





Figure 2: Examples of curves with discontinuous curvature (on the left) and with infinitely smooth 
curvature but still without the Frenet frame (on the right). 



Furthermore, it is a standard hypothesis in the literature about quantum waveguides that the first 
three derivatives of the reference curve T exist and are linearly independent, so that the torsion and 
the distinguished Frenet frame exist. However, this is meaningful only for curves which are three times 
differentiable and have nowhere vanishing (differentiable) curvature n. We find the latter as a very 
restrictive requirement, even for infinitely smooth curves (c/ Figure^]). Indeed, the torsion r is not well 



defined for such curves, so that the limit ([IT]) with the effective Hamiltonian H^ s is meaningless. Partial 
attempts to overcome this technical condition can be found in [3j [5] . In this paper we provide a complete 
answer by considering waveguides built along any twice differentiable curves, with the boundedness of k, 
being the only hypothesis. Our assumptions are very natural and in fact intrinsically necessary for the 
construction of the waveguide as a regular Riemannian manifold. 

Finally, the Gamma-convergence method of [H|6], which seems to work under less restrictive regularity 
once the technical difficulty of the non-existence of the Frenet frame is overcome, implies only (unless 
the waveguide is bounded [S]) a strong-resolvent convergence for ( |1.1| . Furthermore, it does not provide 
any information about the convergence rate. In addition to the regularity issues mentioned above, our 



goal is therefore to use operator methods instead of the Gamma-convergence, establish (1.1) in the 
norm-resolvent sense and get a control on the convergence rate. 



1.2 The content of the paper 

The organization of this paper is as follows. 



In the following Section [2] we explain our strategy to handle the singular limit under mild regularity 



hypotheses and state the main result of this paper (Theorem 2.1). 

We postpone a precise definition of a simultaneously twisted and bent waveguide and of the associated 
Dirichlet Laplacian till Section [3j For reasons mentioned above, we construct the waveguide by using an 
alternative frame defined by parallel transport along the curve instead of the usual Frenet frame. Since it 
seems that this frame is not as well known as the Frenet one, and since we want to include more general 



curves than those usually considered in differential geometry, we decided to include Section 3.1 where 
we thoroughly describe the construction of the frame under our mild regularity conditions. 

The main idea of the present paper consists in smoothing non-differentiable quantities by means of 
the so-called Steklov approximation (see (2.5) below). This procedure is in detail explained in Section El 



The proof of Theorem |2.1| is given in Section [5j Since it is rather long and technically involved, we 
divide the proof into several auxiliary lemmata and the section into corresponding subsections. 
The paper is concluded in Section [6] by discussing optimality of our results. 

2 Our strategy and the main result 

Our strategy how to achieve the objectives sketched in Introduction is based on the following ideas: 

(I) Use the frame defined by the parallel transport instead of the Frenet frame. This alternative frame 
is known to exist for any curve of class C 2 , cf 0. We generalize the construction to the curves 
that merely belong to the Sobolev space Wj ' c °°. 

(II) Work exclusively with the quadratic forms associated with the operators. More specifically, we 
adapt the elegant method of Friedlander and Solomyak [T31 [T5] to deduce the norm-resolvent 
convergence from a convergence of quadratic forms. 

Even if one implements these ideas, the standard operator approach to the thin-cross-section limit 
in quantum waveguides (see, e.g., [5]) still requires certain differentiability of curvature k (which is just 
bounded under our hypotheses). To see it, we sketch the standard strategy now. 

First, one uses curvilinear coordinates, which induce the unitary transform 

Ui : L 2 (Q £ ) -> L 2 (l x uj,e 2 h(s,t)dsdt), (2.1) 

where the Jacobian e 2 h is standardly expressed in terms of k and 6p. In our more general setting enabled 
by the strategy (I) above, we have 

h(-,t) := 1 — et\ (k\ cos 9 — fc 2 sin.9) — ei 2 (k\ sin + fc 2 cos 9) , (2-2) 

where k\ , fc 2 are curvature functions computed with respect to our relatively parallel frame and 9 is 
an angle function defining the rotation of the cross-section euj with respect to this frame. We have 
k 2 = k\ + k\ and, if the Frenet frame exists, our frame is rotated with respect to the Frenet frame by 
the angle given by a primitive of torsion r (cf (|3.4[) below). Consequently, in our more general setting, 



the difference Of — r in ( 1.1 1 is to be replaced by 9 and the effective Hamiltonian reads 

,.2 



tf cff := -A{, - ^ + Oj 2 . (2.3) 



We emphasize that this operator coincides with H* s introduced in (1.1) if T possesses the Frenet frame 



but, contrary to H^ s , it is well defined even if the torsion r does not exist. 



Second, to recover the curvature term in the effective potential of ( 1.1 1, one also performs the unitary 
transform 

U 2 : L 2 (I x u,e 2 h(s,t)dsdt) -> L 2 (I x lu) : U)> H- sVhipl . (2.4) 

The composition U := UiU\ clearly identifies the geometrically complicated Hilbert space L 2 (il 6 ) with 
the simple L 2 (I xio). The standard procedure consists in transforming — A. D S with help of U to a unitarily 
equivalent operator on L 2 (I x w) and prove the norm- resolvent convergence for the transformed operator. 
However, U does not leave the form domain W ' (I x lu) invariant if /ci,fc 2 are not differcntiablc in a 
suitable sense. 

The last difficulty is overcome in this paper by the following trick: 



(Ill) Replace the curvature functions in (2.2) by their e-dcpendcnt mollifications (fi £ {1,2}) 



W-^ .,., *m(0«. ( 2 -5) 



^ r 

8(e) J S _*M 

where 8 is a continuous function such that both 8(e) and e<5(e) _1 tend to zero as e — > 0. 

Then everything works very well (although the overall procedure is technically much more demanding) 
because the longitudinal derivative of the mollified h involves the terms ek^ which vanish as e — > 0, even 



if fcfj diverge in this limit. In more intuitive words, (2.5) can be understood in a sense as that the curve 



is smoothed on a scale small compared to the curvature of the curve, but large compared to the diameter 
of the cross-section of the waveguide. 



The mollification (2.5) is sometimes referred to as the Steklov approximation in Russian literature 
(see, e.g., [I]). At a step of our proof, we shall also need to mollify the derivative of the angle function 9. 

Before stating the main result of the paper, let us now carefully write down all the hypotheses we 
need to derive it, although some of the quantities appearing in the assumptions will be properly defined 
only later. 

Assumption 1. Let T : I — > R 3 be a unit-speed spatial curve, where the interval I C R is finite, 
semi-infinite or infinite, satisfying 

(i) reW 1 2 ' c 00 (7;R 3 ) and k := |f | £ L°°(I). 

Further, let lu be a bounded open connected subset of R 2 and let 9 : I — > R be the angle describing the 
rotation of the waveguide cross-section euj with respect to the relatively parallel adapted frame constructed 
along V satisfying 

(ii) 9 £ W£.°°(J) and 6 £ L°°(I). 

Finally, we assume 

(iii) Vi £ does not overlap itself for all sufficiently small e. 

The conditions stated in Assumption [T] are quite week and in fact very natural for the construction of 
the waveguide f2 e and for obtaining reasonable spectral consequences from ( 1.1 1 (cf Section p^ for further 



discussion). Unfortunately, for making our strategy to work in the case of unbounded waveguides, we 
also need to assume the following (seemingly technical) hypothesis. 

Assumption 2. For any f £ L^ C (I), let us define 



cr f (8(e)) := sup 



\ 



■ n+l 

sup / \f(s)-f(s + V )\ 2 d S . (2.6) 



where e i— > 8(e) is some continuous function vanishing with e. To give a meaning to (2.6) for I ^ R, we 



assume that f is extended from I to R by zero. We make the following two hypotheses 

lim a k (8(e)) := lim £ a^(8(e)) =0, (2.7) 

^=1,2 

lim a 6 (8(e)) =0. (2.8) 

for some positive continuous functions 8, 8 satisfying 

lim 8(e) = , lim -^- = , lim 8(e) = . (2.9) 

e^O e^O 8(e) £^0 

Assumption [2] is satisfied for a wide class of reference curves T and rotation angles 9. First of all, 
let us emphasize that it always holds whenever I is bounded. Indeed, this is a consequence of the more 
general fact that Assumption [2] holds provided that (the extensions of) the representants of / are square- 
integrable functions on R. As other sufficient conditions which guarantee the validity of Assumption [2j 
let us mention that it holds whenever the representants are either Lipschitz, or just uniformly continuous, 
or periodic, etc. In any case, it is a non-void hypothesis for unbounded I only, when it becomes important 
to have a control over the behaviour of ki , ki and 9 at infinity. 

Now we are in a position to state the main result of this paper. 



Theorem 2.1. Let AssumptionU\ and Assumption[Ehold true. Then there exist positive constants £q 
and C such that for all e < Eq, 

U{-A%' - e- 2 E 1 - iy 1 ^ 1 - {H cS - iy 1 @ 1 - 

<c(s + S \\kl\\+e\\k{\\+o- k (5(£))+crM £ ))) > ( 2 - 10 ) 
where O 1 " denotes the zero operator on the orthogonal complement of the span of {if (?) J\\^p £ L 2 (I)} 



and U = U 2 XJ\ is the unitary transform composed of (2.1) and (2.4 1 



Recall that the quantities e||fef||oo an d e||fc|||oo from the right hand side of (2.10) tend to zero as 



0. Hence Theorem 2.1 indeed implies the norm-resolvent convergence of the type ( 1.1 1 and it covers 
all the known results, and much more. Furthermore, the right hand side of (2.10[) explicitly determines 



the decay rate of the convergence ( 1.1 1 as a function of the regularity properties of k\, k 2 and 9. Again, 
it reduces to the well known (see, e.g., |llj ) e-type decay rate for (uniformly) Lipschitz k\,k 2 and 9 
(cf Section ^. 

Assumption ^1 actually requires that the curvatures k\ , k 2 and 9 are not oscillating too quickly at 
infinity (if / is unbounded). We leave as an open problem whether it is possible to have the norm- 
resolvent convergence without this hypothesis. 

We refer to Section [6] for further discussion of the optimality of Theorem |2.1| 

3 Preliminaries 

In the following subsection we introduce the notion of relatively parallel adapted frame for any weakly 
twice differentiable spatial curve. It is then used to define the tube fl e in Section [3T2| while the associated 
Dirichlet Laplacian is eventually introduced in Section [3. 3| 



3.1 The relatively parallel adapted frame 

We closely follow the approach of Bishop [2] who introduced the relatively parallel adapted frame for 
C 2 -smooth curves. Indeed, the extension to curves which are only weakly differentiable requires rather 
minimal modifications. 

Given an open interval / C R (finite, infinite or semi- infinite), let T : / — > M 3 be a C^-smooth 
immersion. Without loss of generality, we assume that the curve T is unit-speed, i.e. |T(s)| = 1 for all 
s £ /. Then T := T represents a continuous tangent vector field of T. 

A moving frame along T is a triplet of differentiable vector fields e^ : / — > M 3 , i = 1, 2, 3, which form 
a local orthonormal basis, i.e., ei(s) ■ Cj(s) — 5ij for all s £ /. We say that a moving frame is adapted 
to the curve if the members of the frame are either tangent of perpendicular to the curve. The Frenet 
frame (if it exists) is the most common example of an adapted frame, however, in this paper the so-called 
relatively parallel adapted frame (RPAF) will be used instead of it (since it always exists). 

We say that a normal vector field M along T is relatively parallel if its derivative is tangential, i.e. 
M x T = 0. Such a field can be indeed understood as moved by parallel transport, since it turns only 
whatever amount is necessary for it to remain normal, so it is a close to being parallel as possible without 
losing normality. 

The RPAF then consists of the unit tangent vector field T and two unit normal relatively parallel and 
mutually orthonormal vector fields Mi,M 2 . Let us note that for any relatively parallel normal vector 
field M, we have (|M| 2 )' = 2M ■ M = 0. At the same time, [M x ■ M 2 )' = 0. That is, the lengths of 
the relatively parallel normal vector fields and the angle between them are preserved. Consequently, the 
definition of RPAF makes sense and it indeed represents an adapted moving frame. 

The existence of RPAF for any C 2 -smooth curve is proved in [2]. However, such a regularity implies 
that the curvature k := \T\ is continuous, which is still a too strong assumption for us. Hence, here we 
provide an extension of the construction of RPAF to curves which are merely T £ IF lo ' c °°(/;E 3 ). This 
implies that the curvature k is locally bounded only, which does not restrict our results whatsoever, since 
the stronger assumption n £ L°° (I) will have to be assumed for other reasons anyway. 



Proposition 3.1 (Existence of RPAF). Let T e W£™(I;M. 3 ) be a unit-speed curve and let Mf and M$ 
be two unit normal vectors at a point T(sq) such that {T(so),M® j-M^} * s an orthonormal basis of the 
tangent space Tr( So )R 3 . Then there exists a unique relatively parallel adapted frame {T, Mi, M2}, such 
that Mi(sq) = Mj 1 and M 2 (sq) = M 2 . The vector fields in this frame are continuous and their weak 
derivatives exist and are locally bounded. 

Proof. By T e W£' C °°(I;R 3 ) we mean precisely that T l £ W£* C °°(I) for i = 1,2,3, which yields that 
r* e C 1 (I) and the derivative T l is locally Lipschitz continuous. This allows us to introduce a continuous 
unit tangent vector field T := T as before and we know that the weak derivative of T exists and is locally 
bounded. 

It remains to find the two relatively parallel normal vector fields M\,M 2 . First of all, let us notice 
that the uniqueness is trivial: the difference of two relatively parallel normal vector fields is also relatively 
parallel, hence preserves the length. So if two such coincide at one point, their difference has constant 
length zero. 

In the first step we find some auxiliary unit normal vector fields satisfying the initial conditions, i.e. 
the vector fields N\,N 2 satisfying N^ ■ T = 0, N^ • N„ = 5^ and A M (s ) = M°, with /u,,v = 1, 2. Such 
fields can be constructed locally by employing the continuity of T and local boundedness of k. Explicitly, 
assuming without loss of generality that one of the coefficients T 1 (so) or T 2 (so) is greater or equal to 1/3, 
we can choose, for instance, 

/ _ T 2 T l \ 

N x := — — , — =,0 , N 2 :=TxNi. 

\V( T1 ) 2 + ( T2 ) 2 v/(Ti) 2 + (T 2 ) 2 ' J 

By means of the fundamental theorem of calculus, we can easily establish the inequality 

\T(s) - T{s )\ <\s- s | ||«||l~((« ,»)) I 3 - 1 ) 

for every s £ J, which shows that iVi, N 2 are well defined in a bounded open interval J around sq. From 
the dependence of the components of JV M on T, we deduce that both iV M £ W 1 '°°{ J; K 3 ). 

In the second step we have to realize that it is always possible to find a continuous and in the 
weak sense differentiable function ■& : J — > K satisfying 'd(so) = and such that the normal vector 
field Mi := iViCOS^ + JV^sini? is relatively parallel in J. This is easily established by expressing the 
derivative of the triple {T, iVi, N 2 } by means of an antisymmetric Cartan matrix and by choosing ■& as 
primitive of the non-tangential coefficient of the matrix coming from the derivatives of Ni,N 2 . Then 
also M 2 := -Nt sini? + N 2 cosz? is relatively parallel, both M M G VF^°°(J;IR 3 ) and Af M (s ) = M°. 

Finally, to get the global existence on /, we can patch together the local RPAFs, which exist in a 



covering by bounded intervals because of (3.1 ). The regularity at the points where they link together is 



a consequence of the uniqueness part. □ 

If {T, Mi, M 2 } is a RPAF, we have 

/ fci k 2 \ ( T \ 

(3.2) 
\ M 2 ) \ -k 2 / \ M 2 J 

Due to r e W '°°(/;R 3 ), functions k\ and k 2 are locally bounded, however they do not need to be 
neither differentiable nor continuous. 

Analogous functions for the Frenet frame, i.e. the curvature k and torsion t, are uniquely determined 
for a non-degenerate curve (i.e. k > 0). Let us examine the uniqueness of fci, k 2 in our general situation. 




Proposition 3.1 says that for a given curve, RPAF is unique if the initial vectors M^, M 2 at some point sq 
are specified. For rotated initial vectors M° := Y^u^i^-v-M®, with \i = 1,2, where 1Z is any constant 
2x2 orthogonal matrix, a different RPAF is obtained in general. Consequently, the functions fc M transfers 
to fc M = Xli/=i ^■/u'^! with /i = 1, 2. Hence the curvatures k\, k 2 are not unique for the curve. 
On the other hand, we have 



K=\T\ = \k1M1 + k 2 M 2 \ = Jk{ + fc 2 , (3.3) 



hence the magnitude of the vector (k\, k 2 ) is independent of the choice of RPAF. Finally, let us assume 
that the curve T possesses the distinguished Frenet frame and let us denote by N the principal normal 
and by B the binormal. It is easy to check that the pair of vectors {Mi, M 2 } is rotated with respect to 
{N, B} by the angle 



0(s) = A) + / r(0 dt 



(3.4) 



where i?o is the angle between vectors Mi(so) and N(so). Consequently, r = fl. Writing, (fcijfe) = 
(k; cos 1?, k sin #) , we can conclude that k and an indefinite integral of r represent polar coordinates for 
the curve (fci, fea), as pointed out in [2]. 

In Figureplwe can see how the pair of Frenet normal vectors {N, B} versus relatively parallel normal 
vectors {Mi, Ma} move along a helix. The longer side of the rectangular cross-section corresponds to 
the direction of N (on the left) and M\ (on the right), whereas the shorter side is the direction of B 
and M 2 , respectively. In the bottom of the Figure, the two frames coincide. 





Figure 3: A waveguide with rectangular cross-section built along a helix. In the left figure the cross- 
section moves according to the Frenet frame, i.e. 9f — 0, so that the waveguide is twisted because 
6 = —t 7^ 0. In the right figure the cross section moves according to the RPAF, i.e. 6 — 0; we say that 
such a waveguide is untwisted. 



Remark 3.2. In Assumption^ of Theorem \2.1\ we state some requirements on the curvature functions 
k\ , k 2 that are not uniquely determined for the reference curve, as we have seen in this subsection. 
However, let us fix some particular RPAF with curvatures k\ , k 2 and recall that the curvatures for different 
RPAFs are only the linear combination of k® and k 2 . When we examine the condition | |2.7| ), we easily 
find that ifki,k% satisfy it, then all their linear combinations do satisfy it as well (due to the triangle 
inequality in L ). Hence there is no ambiguity in Theorem 



2.1 



3.2 The geometry of the tube 

As mentioned in Introduction, the tubes we consider in this paper are obtained by translating and 
rotating a two-dimensional cross-section along a spatial curve T. This definition can be formalized by 
means of the RPAF {T, Mx,M 2 } found in the previous section. 

The cross-section of our tube can be quite arbitrary. We only assume that u is a bounded open 
connected subset of M 2 . The boundedness implies that the quantity 



sup |i| 

tew 



(3.5) 



is finite. We say that u is circular if it is a disc or an annulus centered at the origin of M 2 (with the 
usual convention of identifying open sets which differ on the set of zero capacity) . 
Given an angle function 9 € ^ (/), let us define a rotation matrix 

q _ ( cos 9 sin< 
I — sin 9 cos ( 

Then we define a general moving frame {M®,M%} along T by rotating the RPAF {Mi,M2} by the 
angle 9, i.e., 

2 

M^ = J2K^ M -' A* = 1,2. 

Let f^o := I x cj be a straight tube. We introduce a curved tube $7 e of uniform cross-section ecu as 
the image 

n e :=£(O ), (3-6) 

where the mapping C : flo — » K 3 is defined by 

2 

£(«,t):=r(»)+ e ^t^. (3.7) 

We say that the tube J7 e is &eni if the reference curve T is not a straight line, i.e. n ^ 0. We 
say that O e is untwisted if w is circular or the cross-section is moved along T by a RPAF, i.e. 9 = 0; 
otherwise the tube is said to be twisted (for example of twisted and untwisted tube see Figure pi. A list 
of equivalent conditions for twisting can be found in |17j . 

It is usual in the theory of quantum waveguides to assume the tube £l E is non-self-intersecting, i.e., 
C is injective. The necessary (but not always sufficient) condition for the injectivity is the non-vanishing 
determinant of the metric tensor 

Here <9,; denotes the partial derivative with respect to the i th variable, where the ordered set (s, t\, t 2 ) 
corresponds to (1,2,3). Employing (3.2), it is straightforward to check that the matrix G = (Gij) reads 



where 



We have 



h 2 +e 2 {h 2 2 + hl) -e 2 h 3 -e 2 h 2 
G = | -e 2 h 3 e 2 | , (3.8) 

-e 2 h 2 e 2 



h(-, t) := 1 — et\ (fei cos6> + k 2 smQ) — et 2 {—k\ sin# + k 2 cos 9) , 
h2(;t):=-tj, (3.9) 

h 3 (;t) :=t 2 e. 

\G\ :=det(G)=e 4 /i 2 , (3.10) 



hence the condition on the determinant being everywhere nonzero requires that h is a positive function. 
The latter can be satisfied only if the functions hi , k 2 are bounded. Therefore we always assume 

K€L°°(I), (3.11) 



which is equivalent to the boundedness of fei, k 2 due to (3.3). In particular, we have 

HMoo < Moo < oo, /*= 1,2, (3.12) 

where || • H^ := || • ||j,<»(7). Using in addition the boundedness of w, we find the bound 

h(s,t) > 1 -eoll/clloo (3.13) 

for every (s,t) £ fio- This ensures the positivity of h for all sufficiently small e. 
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If the determinant (3.101 is positive, the matrix (3.8) is invertible and we have 



G~ 




h 3 

~ 2 h 2 + 



h\ 




Summing up, assuming (3.11 1 and the injectivity of C, the mapping induces a global diffcomorphism 
between the straight tube f^o and il e , and the latter has the usual meaning of a non-self- intersecting 
curved tube embedded in M 3 . For sufficient conditions ensuring the injectivity of C we refer to [TO) 
App. A]. The above construction gives rise to Assumption [11 

Remark 3.3. Abandoning the geometrical interpretation of Q s being a non-self-intersecting tube in M 3 , 
it is possible to consider (f2o,G) as an abstract Riemannian manifold, not necessarily embedded in M 3 . 
This makes ( |3.11| ) (together with the smallness of e to ensure that the right hand side of (3.13) is positive) 
the only important hypothesis in the present study. In other words, the injectivity assumption (Hi) in 
Assumption^ can be relaxed, the results of the present paper hold in this more general situation. 



3.3 The Hamiltonian 

Let us now consider Q e as the configuration space of a quantum waveguide. We assume that the motion 
of a quantum particle inside the waveguide is effectively free and that the particle wavefunction is 
suppressed on the boundary of the tube. Hence, setting h = 2m = 1, the one-particle Hamiltonian acts 
as the Laplacian on L 2 (Q, £ ) subjected to Dirichlet boundary conditions on dil e : 

-Ag E . (3.14) 

The objective of this subsection is to give a precise meaning to this operator. 

The most straightforward way is to assume that C is injective and define (3.14) as the Dirichlet 
Laplacian on L 2 (il e ). Indeed, this is well defined for open sets and from the previous subsection we 
know that C induces a global diffeomorphism, so that, in particular, f2 e is open. More specifically, the 
Dirichlet Laplacian — A^ is introduced as the self-adjoint operator associated on L 2 {p, e ) with the closed 
quadratic form 



Q£M:=l|V^||i a(n , 



Dom(Q^) = <' 2 (fi £ ). 



From this point of view, we regard the tube as a submanifold of R 3 . For the description of J7 e , the 



most suitable coordinates are the curvilinear 'coordinates' (s,t) € f2 defined via the mapping C in (3.7). 
They are implemented by means of the unitary transform 



Ut : L 2 (n € )^H e := L 2 (n Q ,\G(s,t)\ 1/2 dsdt) : {if} i-> %j) o £} 



(3.15) 



mentioned already in (2.1). The transformed operator H e := Ui(—A D E )U 1 can be determined as the 



operator associated with the transformed form 



1,2/ 



Dom(Q £ ):=£W ( T(ft £ ) 



(3.16) 



Here G %3 are coefficients of the inverse metric (3.8) and the Einstein summation convention is adopted 



(the range of indices i,j being 1,2,3). In a weak sense, H £ acts as the Laplace-Beltrami operator 
— |G| -1 / 2 9j|G| 1 / 2 G' J '9j, but we shall not need this fact, working exclusively with quadratic forms in this 
paper. 

Let us emphasize that, for the quadratic form Q £ to be well defined, the matrix G does not need to 
be differentiable, a local boundedness of its elements is sufficient. As a matter of fact, the form domain 
U\W ' (ri ) can be alternatively characterized as the completion of G^°(f2o) with respect to the norm 

U\\q € ■= yJidi^&id^ + M^. 

If the functions k and 6 are bounded, it is possible to check that the Q e -norm is equivalent to the usual 



norm in W ' (£7q)- For this reason, in addition to (3.11), we assume henceforth the global boundedness 



6e L°°(I). 



Then we have 



Dom(Q e ) 



w 1,2 (no) 



(3.17) 
(3.18) 
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Remark 3.4. Now, if C is not injective, the image (3.6) might be quite complex and the standard 



notion of the Dirichlet Laplacian on Q e meaningless. Nevertheless, the Laplace- Beltrami operator on the 
Riemannian manifold (f?o, G), i.e. the operator associated on T~i e with the closure of the form defined via 
the second identity in (3.16) on the initial domain C^(flo), * s fully meaningful. Moreover, it coincides 
with H e defined above. This is the way how to transfer the results of the present paper to the more 



general situation of Remark 3.3 In particular, Theorem 2.1 holds without the injectivity assumption (Hi) 
in AssumptionUj provided that we properly reinterpret the meaning of (3.14) as the Laplace- Beltrami 
operator H e in (fioi G) and we write just TJ<i instead of U in the statement of the theorem when dealing 
with the more general situation. 

Finally, let us recall that the spectrum of H £ explodes as e~ 2 E\ in the limit as e — > 0. This is 
related to the fact that the ground-state eigenvalue of the cross-sectional Laplacian — A £ p equals e~ 2 Ei. 
Therefore, to get a non-trivial limit, we rather consider the renormalized operator 



H, 



H F 



i E 1 



in the sequel. 



4 The mollification strategy 

Our strategy how to reduce the regularity assumptions about the waveguide consists of the three points 
(I)— (III) roughly mentioned in Section [2| T he first of them, i.e. the usage of RPAF instead of the Frenet 
frame, was already explained in Section 3.1 The item (II) consists in working with associated sesquilinear 
forms instead of operators. In the preceding Section |3.3| we introduced the Dirichlet Laplacian in the 
tube Q e using exclusively quadratic forms and it enables us to understand the derivatives in the weak 
sense and to reduce the requirements on the differentiability of the reference curve. 

However, as explained in Section [2] for the standard operator procedure to work, certain additional 
smoothness of curvature functions are still needed. In this section we propose a method how to proceed 
without any extra regularity hypotheses. It is based on the mollification procedure (III) sketched in 
Section [2] and we believe it might be useful in other problems as well. 



4.1 The modified unitary transform 

As explained in Section [2j the main idea consists in mollifying the curvatures k\ , ki by means of the 
Steklov approximation to get fef , fc| introduced in (2.5). If / is finite or semi-infinite, we adopt the 
convention of Assumption [2] to give a meaning to function values outside / in the definition. That is, we 
assume that fc M , with fi = 1,2, are extended from / to R by zero. 



The definition (2.5) involves a positive continuous function e *— > 5(e) which is supposed to satisfy 

and lira -^ = 0. (4.1) 



lim 5(e) 

£->0 



o 5(e) 



Here the first assumption is reasonable since then fcfj > fc M in a certain sense (see Section 

The relevance of the second condition will become clear in our computations. 



4.2 



It follows from (3.3) that 



1,2. 



Furthermore, the mollified functions are different iable for any positive e, 



K(») 



kfi(a- 



He} 
2 



-kp(s- 



S(e) 
2 



5(e) 



M=l,2, 



below). 

(4.2) 

(4.3) 



although the derivative might diverge in the limit as e — > for non-differentiable functions k^ 



We also introduce a smoothed version of the 'Jacobian' (3.9) 

h e (-,t) :=l-sti (kl cos 8 + k^sm 8) - et 2 (-fcf sin6> + k\ cosO) 
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and of the determinant ( 3.10[ ), \G\ := e 4 h 2 . The latter will be used to generate a modified version of the 
standard unitary transform U% from Section [2] We define 



u 2 ■ n e -> n e ~ l 2 n 



\G(a,t)\V* 



ds dt 



{^^IGI^V} 



(4.4) 



The norm and inner product in the Hilbert space H £ will be denoted by 



and (',') e , respectively. 



In addition to (3.111, let us assume (3.17) in the following, so that the form domain of H is given 
by (3.18). Since the Sobolev space W ' (f2o) is left invariant by the modified transform U 2 , the operator 



H s := U 2 H s U 2 l = \G\ 1/4 H e \G\- 1/4 
is well defined in the form sense. The associated quadratic form reads 



QebP] 



1 

n hh 
1 



1 



ds dt 



1 



e 2 ./On h e 



'V'ijj\ 2 dsdt 



£ 2 ./n h e 



\ip\ 2 dsdt 



(4.5) 
(4.6) 



., (k 1 k s 1 + k 2 k e 2 M 2 dsdt-- 

^ -'Ho e 



n h e 



{{k\f + {klf)\^dsdt 



(d s + 0d a )h t 



Ahhl 



-\tp\ 2 ds dt - 



(d a + 6d a )h E 

hh 2 



Re(ip(d a + 9d a )ip)dsdt, 



-1,2 



with tp G Dom(Q £ ) = W Q ' (CIq)- Here V is the gradient operator in the 'transverse' variables (£1,^2) 
and d a is the transverse angular-derivative operator 

/ n „; d d 

d a := (t 2 , -h) • V = £ 2 — - t 17 — . 
oti ot 2 

An important feature of the operator H e is its boundcdncss from below. We prove it together with 
another relation used in our computations below. 



Lemma 4.1. Let Q e be the quadratic form defined by (4.6) and let (4.1) be satisfied. Then for all 
ip G W ' (SIq) and small enough e 



QM> -2J a ,k h . 

Proof. If we assume that e is so small that 



(d s + 6d a )ip ds*- 91^0 VII 2 . 



(4.7) 



3 5 

- < 1 -eallKJIoo < h < 1 + ea||K|| 00 < - , 



then the same relation holds for h e and using (3.12) we easily get 



\f ^(kiki + k 2 k^\ 2 dsdt-lj^ ~{{kt? + {ktY)^\ 2 dsdt ^sMUnl- 

The estimate on terms proportional to e~ 2 is based on the Poincare-type inequality 

I \V(p\ 2 dt~E 1 I \<j)\ 2 dt>0 

that holds for all <f> G W ' (w), and on Fubini's theorem. Due to nontrivial measure in our integrals, we 
have to use substitution <\> :— y/h/h E ip, then we obtain 



1/ A|vvi 2 d S *-§ / htfdsdt 

£ Ja n e £ Jn n e 

= - 2 I (|V^| 2 -Si|^| 2 ) dsdt 
£ Ja 



3(fc 2 + k 2 ) kikf + k 2 k e 2 (fcf ) 2 + (fef) 



■„£\2 



>-3||«| 



4/i 2 2/i/i E 4/i 2 

|0| 2 d s dt = -3|| K || 2 o ||V|| £ 2 . 



|</>| 2 dsdi 
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The last term in (4.6 ) can be estimated using the Schwarz inequality and the simple Young's inequality 
(2ab <a 2 +b 2 ) 



(d 3 + 9d a )h £ 
hh 2 



Re(ip{d s + Od a )ip) ds dt 



< 



(d. + 9d a )h t 
Ahhi 



-\ip\ 2 dsdt- 



hh F 



(d s + 9d a )i; 



ds dt 



Summing up, we get 



{d s + 6d a )?P dsdt -8||/s|£,|M|; 



(d s +9d a )h e y 

4hM 



\i/}\ 2 dsdt 



and the proof is finished by the estimate 

2 



((d s + dd a )h ( 

4hM 



-\ip\ 2 dsdt 



<51 



""' a 2 \\4lW\ 2 s<\\4loU\\ 2 e 



8(e) 



which holds for small enough e due to (4.1) 



□ 



As a consequence of Lemma 4.1 we get that for a constant A < — 9||k||oo the operator H £ — A is 
invertible and it holds 

^- 1 " - ' (4.8) 



\(H e - A) \\ B (H e ) < 



|A|-9|| K ||L' 



4.2 Convergence properties of the Steklov approximation 

Let / be a bounded function defined on an interval I and let f e be its Steklov approximation 

i r + ^ 



(4.9) 



where S is a positive continuous function on R satisfying the first of the requirements of (4.1 1. Again, 
we recall the extension convention of Assumption [2] if / is finite or semi-infinite. 

In the computations below we require for / = ki,k 2 ,0 that in a certain sense f e converges to / in 
the limit e — > 0. Namely, the crucial requirement reads 



Ju-r\ 2 \v\ 2 d S 



1/2 



^>0, W>e<< 2 (/). 



(4.10) 



In the following we will find the requirements on / such that this condition is satisfied. 



We shall start with the estimate on the left hand side of (4.10) 



Lemma 4.2. Let f € L°°(I) and let f £ be the Steklov approximation of f. Let ip <E W l,2 (L) and 
finally let {a n } n = n _ C I . n G Z be the strictly increasing sequence of numbers where a n _ — inf se / s, 
a n+ — sup se/ s, all the intervals I n := (a n , a n+ \) are finite and n± can be either finite number or ±oo. 
Then 



\f - r\ 2 \<p\ 2 ds < sup 

/ ?l_<n<Tla 



|/-/ £ " 2 



L 2 (/J 



■2||/-/ E||:j 



L 2 (I„) 



O-n+1 — 0"n 

Proof. The main idea of the proof is rewriting the estimated expression as 



m2 



(4.11) 



[\f~f\ 2 M 2 d S = J2 r" 1 9eM 2 ds= J2 {[9 n eW\XT - r fl 5?2Re(w)d S 
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where we integrated by parts and where we defined, for all s G /, 



Xn{C) i s the characteristic function of the interval I n (i.e., Xn(£) = 1 f° r £ € [onjGn+i] and Xn(0 = 
elsewhere) . The proof is then completed by using the Schwarz and Young inequalities and the relations 



sEl 



a„+i 



supg^s) = g?(a n+1 ) = \ |/(£) - f(i)fdi = ||/ - f 6 \\h (In) , 



\<p(a n+1 )\ 2 < 



0>n+l 



+ 1 \\Lp\\ W l,l(I n ). 



D 



The lemma is of great importance for our computations. From the generalized Minkowski inequality 
it follows that 



11/- / £ ||l2(/„) 



< 



sup 



kl<- 



1/2 

\f( s )-f( s + r j)\Us) =:u 2 (8(e), /,/„). 



(4.12) 



Here the notations! symbol lo p is adopted from [T], where it is referred to as modulus of continuity 
generalized to space LP and is computed for a function /, positive number 8(e) and interval I n . In pQ it 
is also shown that this quantity tends to zero with 8(e) if the interval I„ is finite, rl This directly yields 
that if the interval I is finite, the convergence (4.10) holds. 



If I is infinite, we can cut it into finite intervals (a n ,a n+ i) and for each n the expression in square 



brackets on the right hand side of (4.11 1 would tend to zero when e — ► 0. Unfortunately, the supremum 



over n_ < n < n + might not have the zero limit. On the other hand, this may happen only in a case of 



functions that oscillate quickly at infinity (see Example 6.2). In other words, also for unbounded i", the 



condition (4.10) is satisfied for functions that behave 'reasonably' at infinity. 



We summarize the above ideas in the following proposition. For simplicity, we use the result of 
Lemma |4~2| for the equidistant division a n := n and we recall the extension convention of Assumption [2] 
for finite or semi-infinite /. 

Proposition 4.3. Let ip <E W ' (I), f E L°°{L) and let f e be the Steklov approximation of f given 
by (pL9|. Then 

1/2 



b'- r 



lipids 



< 



X 3 o- f (5(e)) \\ip\\ w i,2 {1) 



(4.13) 



where o~f(8(e)) is given by (2.6). Furthermore, the quantity a f (6(e)) tends to zero as e — ¥ if any of the 
following conditions is satisfied: 

(i) / is finite, 

(ii) / periodic on L eyA , 

(iii) / G L 2 (L cxt ), 

(iv) / G C°(L^T t ), 

where L cxt '■= I \ L m t for some finite open (possibly empty) interval L mt . 



convergence properties of (2.6) 



Proof. The inequality (4.13) follows from Lemma 4.2 and the relation (4.12). It reminds to prove the 



Any bounded interval / can be covered by intervals /„ = (n,n + 1], n G I with X C Z finite. Then 
a f (8(e)) is proportional to maximum of u>2(8(e), /, /„) over n G I. Since u>2(8(e), f, L n ) converges to zero 



1 More precisely, the proof in [T] is made for u>2(S(e), /, M) =: u>2 (S(s), /) where / £ L 2 (R). Here we consider a bounded 
function defined on finite I n and prolonged onI\/„ in such a way that the extension is an L 2 -function. Then the proof 
from [l] can be used. 
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for every n £ I as we explained above, a f (8(e)) converges to zero if interval / is finite. For the same 
reason the convergence of a f (8(e)) does not depend on the behaviour of / on a bounded subinterval 
Tint C I in case of infinite or semi-infinite I . 

Also the convergence of a f (8(e)) for periodic functions follows from the fact that the period of length 
q < oo can be covered by finite number of intervals (n, n+1]. (More straightforwardly, using the sequence 



a n = nq in Lemma 4.2 we get a f (8(e)) ex uj 2 (8(e), /, (0, q)), which tends to zero as e — > 0.) 

In the case of Z-Munctions, we use the fact that for any e > there exists a finite interval I[ n t such 

that f(s) < e for all s £ M. \ I{ n t. Then again the significant contribution to a f (8(e)) comes from the 

finite interval. 

Finally, for the uniformly continuous functions the situation is even more simple. Here already the 

quantity 

Woo(<J(e),/):= sup |/(6) -/(6)l 6,6e/ ext , 

|&-&|<*(e) 

called the modulus of continuity in [lj, tends to zero when 8(e) tends to zero and a f (8(e)) can be 
estimated by u> 00 (8(e), /). Let us note that if I is bounded or semi-bounded, the function need not to be 
uniformly continuous after extension by zero outside /. However, the convergence can be still ensured 
by dividing I on the (5(£)-neighbourhood of the end point(s) and the rest of /. Then we can estimate 
the integral over inner part by uj oc (8(e), f) since here the function is indeed uniformly continuous, the 
remaining integral can be estimated by ||/||oo<5(£) which tends to zero as well. □ 

5 The norm-resolvent convergence 

In this long and technically demanding section we give a proof of Theorem |2.1| 



5.1 Comparing operators acting on different Hilbert spaces 

We start with describing a way how to understand the resolvent convergence of operators — A D E and H c g 
acting on different Hilbert spaces L 2 (Q, e ) and L 2 (I), respectively. 

First of all, we recall that, in Section |4~Tj we introduced the operator H e on H E which is unitarily 



equivalent to -A D £ — e 2 £' 1 . It is therefore enough to explain the resolvent convergence of H £ and H e g. 
Our strategy is to reconsider these operators as certain operators on the fixed Hilbert space 

n :=L 2 (Q ) 

and to show that the error due to the replacement becomes negligible in the limit as e — > 0. 

Recall that we have denoted the norm and inner product in the e-dependent Hilbert space T-L £ by 
|| • || e and (v) e , respectively. We simply write || • || and (•, •) for the norm and inner product in Hq. 
Finally, || • ||/ and (•, •)/ stand for the norm and inner product in L 2 (T). 

In order to have a way to compare operators acting on H s and Ho, let us introduce yet another 
unitary transform 

IGI 1 ^ 



U e :H s -^n :^^^—i>j. (5.1) 

For the convenience of the reader, we present here the following diagram explaining the relation with 
the other unitary transforms introduced so far: 



L 2 (n e ) Ul > U e u e (5.2) 




It is important to emphasize that while the transformed resolvent U e (H e — i) U s 1 on T-Lo is well 
defined (as a unitary transform of a bounded operator), the similar expression for the (unbounded) 

16 



operator H e may not have any sense. Indeed, H £ acts as a differential operator, while \G\ may not be 
diffcrcntiablc under our minimal assumption. The same remark applies to H e and U2, as already pointed 
out in Section [21 

Summing up, using the unitary transforms described above, it is possible to reconsider the resolvent 
of the Dirichlet Laplacian —A D E as an operator on Ho- It remains to explain how to reconsider H e s 
acting on L (I) as an operator on the 'larger' space Ho- This is done by introducing the following 
subspace of Ho: 

Hl:={^£H I 3ipeL 2 (I), il>{s,t) = ip{s) Ji(i)}. (5.3) 

Recall that J\ denotes the eigenfunction corresponding to first eigenvalue of the transverse Dirichlet 
Laplacian — A^,; we choose it to be positive and normalized to one in L 2 (u>). Hi is closed, hence 

H =Hl® (Hi) 1 - (5.4) 

and every function if) £ Ho can be uniquely written as 

^> = PiV + (1 - Pity =: faJi+^ (5-5) 

with ij)\Ji £ Hi, y ^ {Ho) 1 " and Pi being projection on H\, 

(Piif>)(s,t) := (J Ji(t)if>(a,t)dt\ Ji(t) = Ms)Ji(t)- (5.6) 

To shorten the notation, we denote by ipi^Ti the function ipi®^Ji, i.e. the function on/xw which assumes 
values as ipi(s)J\(t). Such a decomposition of functions ip £ Ho will be extensively used throughout the 
text with the same notation. 

Now we can introduce the isometric isomorphism 

7T : H\ -> L 2 (I) : {Vi (s)Ji(i) >-> lM*)} ■ 
Let q e s be the quadratic form associated with the operator H c g, i.e., 

q cS {<p} = J \^(s)\ 2 ds + C U J 9( S ) 2 \^( S )\ 2 ds -~J ^sf^fds , Dom (q cS ) = W^(I) . 

(Recall that the basic Assumption 111 requires that both k and 9 are bounded functions, so that H e ff is 
well defined as a bounded perturbation of the one-dimensional Dirichlet Laplacian — A^,.) The form q e s 
can be identified with the quadratic form Q c h acting on the subspace H\ as 

QesbPiJi] ■= f \d s ^iJi\ 2 dsdt-- ( K 2 \i>iJi\ 2 dsdt + C u f 9 2 \faJ 1 \ 2 dsdt = q cS [fa], 

Jn 4 Jn Jn (5.7) 

Dom (Q cS ) := {^ e H\ \ fa £ W^ 2 (I)} . 

In a similar way we can identify operators acting on Hi C Ho and L 2 (I). We tacitly employ the 
identification, without writing down the identification mapping n explicitly in the formulae. In particular, 
denoting by 1 " the zero operator on (Hi)- 1 , the operator (H e g — i)~ @ 1 - can be understood as an 
operator acting on the whole space Ho- 

5.2 Proof of Theorem EO 



At first let us explain the connection between the operator U(— A D E — e~ E±— i)~ U from formula (2.10) 
and the operator U e (H e — i) _1 C/ £ _1 we spoke about in the previous section. More precisely, we show that 
these two operators are identical. Indeed, recall that U — U2U1 where U\ and U2 arc unitary transforms 
described in Section^and in addition that H £ = U2Ui(—A D ' — e~ 2 Ei)Ui U^ . Using the diagram (5.2) 
we easily get that 

U 2 Ui(-A$ - £- 2 E x - i)- l U^ l U 2 - 1 = U e U 2 Ui(-A% - e- 2 E l - i^U^U^U- 1 = U e (H e - i^U' 1 



and in following we will prove Theorem 2.1 using the last expression. 
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Another point is that according to [TB] (Theorem IV. 2. 25), if the formula (2.10) is satisfied for a A 
from the resolvent set of H Q g, then it holds true for all such A. In particular there exists a constant C\ 
such that 



U e (H e - i^U' 1 - ([H eS - i) l ®0^ 



8(Ho) 



u e (H E - xy'u- 1 - ((tf eff - A)" 1 



®0 J 



B(tto) 



(5.8) 



Hence our aim is to prove that the right hand side of the last expression tends to zero for some A < 
— 9||t||^o, since such A belongs to resolvent set of H c g and also of H e (cf (4.7 ) which yields Q E [^] > 
-SJkOV'llf, similarly q eS [ip] > - 



2 
I oo I 



/. r^llV'll/). This proof is divided into proof of two auxiliary 



lemmata, where in every lemma we compare one of the operators on the right hand side of (5.8) with 
the resolvent of operator 



Ho 



1 






K 

T 



a 



i. 



(5.9) 



The crucial step lies in comparison of Hq and H e stated in the first of these lemmata. 

9||k;|| 2 be a real constant and let the assumptions of Theorem 2.1 be satisfied. 



Lemma 5.1. Let A < 

Then 



|[/ e (H £ -A)- 1 C/- 1 -( J ff -A)- 



\B(Ho) 



<c 



(e + e(\\k{\\ c 



v2||oo 



a(5(e)) 



(5.10) 



for some constant C . 



Second lemma giving the comparison of Hq and H c s represents only a tiny improvement of the result 
above. 



Lemma 5.2. Let H be the operator defined by (5.9) and let H c g be the effective Hamiltonian (2.3). 
Then 



(H a - xy 1 - {(H cS - xy 1 © o- 1 



\B(U 



<Ce 



for some real constants C and X < — 9||/c| 



-112 
oo * 



Proofs of Lemmata |5.1| and |5.2| will be given in Sections 5.4 and 5.3 respectively. These lemmata will 
be proved using a trick employed originally in |13j , where the estimate on the norm of the difference of 
resolvents is obtained by a usage of the associated quadratic forms. We state the following proposition 
to give the reader basic idea of this trick; in our proofs we have to modify it due to distinctness of Hilbert 
spaces our operators act on and the main idea could be hidden by loads of technicalities. 

Proposition 5.3. Let A, B be positive self-adjoint operators acting on Hilbert space H. with A^ 1 , B^ 1 G 
B(H) and let Qa, Qb be associated sesquilinear forms with DomQ^ = DomQs. Let us assume that for 
all <f), ip G Dom Qa 



\q a (<i>,ii>) - Q B (<i>,n < vVQaMWQbM 



(5.11) 



Then 



IA- 1 -B- 



\B(U) 



<a 



\A- 



\B{H) 



\B- 



\B(H)- 



Proof. Due to the assumption (|5.1l|) it holds that for all /, g G H 

\Qb 



| (/, (A- 1 - B- l )g) | = | (A4>, (A- 1 - B-^Byj) \ = \Q B (<f>, $) - Q A {cf>, rp)\ < 

< <j^qmVQb~W] < ey/U-mBmy/WB-iWBMWfWhW 



(5.12) 



where the choice / = A<f>, g = Bt/j is possible for all /, g G "H due to boundedness of A^ 1 and B^ 1 . 
This choice ensures that <f> G Dom Q^ n Dom A which due to the representation theorem (see [H]) yields 
(A(f>,ip) = Qa(<I>,iP) for all ip G DomQ^. Similarly ip G Dom Qb H Domi? and ((f), Bip) = Qb(4>, VO f° r 
all 4> G DomQg. Inequality (5.12) yields directly the statement of the Proposition. □ 
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5.3 Proof of Lemma 15.21 

The quadratic form associated with the operator Hq reads 



Qo[il>]= f \d s ^\ 2 dsdt + \ f {V^dsdt- 1 ^- [ 



\ijj\ 2 dsdt 



(5.13) 



•a 



' dsdt 



n 2 \ip\ 2 dsdt 



with 



rl.2. 



DomQ = W^(fl ). 
Due to the equality — A£,Ji = EiJi, this form acts on the set W Q ' (fig) H Wj in the same way as Qoff 



rl,2, 



given by (5.7) which we identify with the quadratic form q c g acting on W ' (I). This yields 

qes(<pi,ipi) - Qo(<f>iJi,i>iJi) = 0, 



(5.14) 



where we use the notation from (5.5 1 and which leads us to slightly modified choice of functions f,g 

We assign Hq 



comparing to the proof of Proposition 



5.3 



H 



^> H eft 



H, 



eft 



A and we choose 



/ := Hq4>, Pig := (H^ipi)Ji, with g A -,i))- L unspecified. If we denote by Qq and q* s , the quadratic forms 



associated to Hq and H* s , respectively, we can rewrite the term analogous to the one estimated in (5.12) 



a.s 



(ifoV, [(ftf) l - ((tf c A ff ) '©O 1 )] ((H^)^ + g 1 - 



Qcft( 



l, Vi) + (</>"-, (H^i)Ji) + (<t>,9 X ) ~ QfoiJuipiJi) ~ Qo (<t> X ,i>iJi) = (^ff x ) 



Here all the terms except of (4>,g ± ) vanish due to (5.14) or due to the orthogonality of ip , (p 1 - and J\. 
Hence we can estimate 

(/, [{ElY 1 ((H^y 1 ©0^)] g)\ = 1(0,^)1 = |((ffo A ) _1 /,(l - Pi)5)| 

<ll/IINIII(^r 1 (i-i > i)lk«o) 

< £ C J . v /||(^)-i||||/|||| fl ||, (5.15) 

where the last estimate follows from the relation 

||^||<eC ±v /Q^] 



that will be proved in Section 5.5 The proof is completed using the estimate 

1 



( H o) lle(Wo) 



< 



IAI 



(5.16) 



analogous to (4.8), i.e. with C := C± 



|A|- 



5.4 Proof of Lemma 15.11 



In thi s pr oof the crucial and most tedious point is to check that the assumption (5.11) of the Propo- 
holds true. Let Qg[^] := Qe[ip] — ^IIV'IU be the quadratic form associated with the operator 



sition 
H R - 



5.3 



A =: H* . Recall that Qo was introduced in previous section. If we assume that these forms 
(understood as sesquilinear forms) satisfy 



|Q £ A (0,^)-g^(0,V)|<^(e)v / Qo['/']a A M 



(5.17) 



rl,2 



for all (/>, ip <G W ' (fio); then we can derive the statement of Lemma 5.1 using similar ideas as in 



Proposition |5.3| We only have to realize that the operators H e and Hq act on different Hilbcrt spaces, in 
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fact H e is not self-adjoint on T-Lq where Hq acts. However, the Hilbert spaces T-L e and Ho can be identified 
via the unitary transform U B defined in (|5.1[), which leads to the estimate 



\(f,[U e (H^U- l -(H*)-']g)\ < ^((ftf)- 1 /,^, 



A\-l, 



Q*((H*)- l f,(H^g 



+ | (/, (U e - l^H^U^g) | + | (/, (H^)- 1 ^ 1 - l)g) | + | ((H>)- V, (Of - l).g) | • 
The operator U £ differs from identity only by amount proportional to e, which we use in the estimate of 



last 3 terms. Together with (5.17) the final estimate reads 



|(/, [UeiH^U- 1 - (If*)- 1 ] 5 )| < ( C £ + 2cr( e ) v /||(^)-i|| B(Wo) ||(^)-i|| B(We) ) ll/H ||. 9 || (5.18) 

where c = 12a||Ac||oo (||(-ff^) _1 ||0(« e ) + IK-^o ) _1 |Ib('Ho)) and where we can in addition estimate the norms 

IK^)- 1 !!^), ||(ffo A )"Uk«o) by © and (OS). 

It remains to prove (5.171, i.e. to find a(e) and to prove that this quantity tends to zero provided the 
assumptions of Theorem |2.1| are satisfied. This part of the proof is very technical and lengthy, on the 
other hand, the reasons for our rather complicated assumptions will be explained. 



5.5 Proof of relation (5.17) 



rl,2. 



Let 4>,ip & W : (f2o). We have to establish suitable estimates on the difference of the following sesquilin- 
ear forms 



>.!■>=/ r£-(d. + 9d a )$(d B + dd a )il>d8dt+ 

1 f h _,-_,,, , Ei f h 



(5.19) 



, , —V'<f>-Vihdsdt- -4 
£ 2 ./fin h e e 



-<fnj} ds dt 



On '% 



1 - 3 

—77(k 1 kl + k2k2)4>ipdsdt-- . . 
o /le 4 J no ^ 

(d 8 + 9d a )h e 



({klf + (k e 2 ) 2 )<j)^dsdt 



2hh\ 

(d s + ed a )h t 

AhM 



(c/)(d s + 9d a )ip + (d s + 6d a )<H)) ds dt 



- cfiip ds dt — A / — (fnp ds dt 
Jn h £ 



and 



?o (4>, ip)= d s <j)d s ip ds dt + C u / r# ds dt 

+ ~~2 / V '^ ' V '^ dlS * 1 / ^ ds dt 

£ JQn £ -'On 



1 



(k\ + kl)W> dsdt- A / </V ds dt 



(5.20) 



On the right hand side of the estimates the term \/Qg [4>]Q^ [ip] should stand. However, due to Lemma 
and similar statement on Qg[</>] we get the inequalities 



4.1 



Q^} > \\\(d s + 0d a W + \(-X~ 9||«||L) IIV^II 2 , 

QM>\\dM 2 + (-\- l ^)U\\ 2 . 



(5.21) 
(5.22) 



Recall that we assume A < — 9|I k II?oj hence in front of ||"0|| 2 j ||0|| 2 there stand positive numbers and we can 
come from estimates by \/Qq [0] , \/Q £ W\ to estimates by norms HV'II 2 , ||<^|| 2 or norms like ||(9 s + 09 a )^||, 
ll^ll- " 
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Some of the estimates can be performed easily using just the Schwarz inequality and straightforward 
estimates: 



((d a + 9d a )h l 

AhM 



-(pipdsdt < C X E- 



(ll*fl 



W®)NA> 



(g, + ed a )h e 

2hh 2 



(<p(d s + 9d a )ip + (d s + 6d a )(pip) ds dt 



<Coe 



(\\K\\oo + \\k e 2 \\oo)yfQMQM, 



(5.23) 



A / -—(j/ipdsdt — X I <p?p ds dt 

o he Jflo 



<C 3 eJQ^[4>}QM 



where the constants read C\ 



,C 2 



and C3 



12aA 



The hrst two inequalities are stated 



in terms of the quantities ||fcf ||oo, which can be replaced by ||fc^||oo for sufficiently regular curves T, so 
that the results of previous papers are recovered. On the other hand, for non-diffcrcntiable curvatures 



the Steklov approximations (recall (4.3 1 ) yields 

e (ll*?IU + ||*S 



< 4 



5(e) 1 



where the right hand side tends to zero due to the second assumption in (2.9 1. Summing up, all the 



terms estimated in (5.231 tend to zero 



To estimate the rest of terms on the left hand side of (5.171, the Hilbert space decomposition (5.4) 



has to be used. The following computations will also show why the assumption (2.7) is needed 



5.5.1 Hilbert space decomposition and estimates by a k (5(e)) 

Let us estimate the difference of terms on the third lines of (5.19[) and (5.20). It is easy to find 



k(0»l 



1 - 3 

^(kik 6 + k 2 k 2 ) (pipdsdt - - . 
n h l ~ 4 J Qo hj 



((k e ) 2 + (k £ 2 ) 2 ) (pipdsdt 



1 



(k\ + k\)(t>ip ds dt 



< 



3|NU I (\k 1 -k\\ + \k^-kl\)\ ( pmdsdt + sa\\K\\l U\\U\\ 
Jn 



(5.24) 



Then the first term on the last line can be estimated by the Schwarz inequality to get a product of 

/ \l/2 

the term (J n (|fci — fcf | + \k 2 — fc||) \4>\ 2 ds dt ) and the analogous one with ip instead of <p. However, 
to proceed further, we have to use the Hilbert space decomposition. If we rewrite the function (p (and 



similarly tp in the other term) as in formula (5.5), we get 



f (Ih-kH + lh-k^l^dsdt 

Jq 



(\h - k\\ + \k 2 - k £ 2 \) (|0i| 2 Ji 2 + 2Re {4>iJi4> 



i_L|2\ 



ds dt 



and we will be able to prove that the last integral tends to zero. The term containing |0i| 2 is (after 



another estimate by the Schwarz inequality and recalling the normalization of J\) analogous to (4.10) 



and tends to zero according to Proposition |4. 3| and assumptions of Thcorcm |2.1| The mixed terms vanish 
due to orthogonality of (p and J\ and thanks to the fact that neither fc M nor k 6 depend on the variable t. 
The remaining part with \<p \ 2 tends to zero according to the following ideas. 
Using straightforward estimates it is possible to find 



a,,,. 1 1 - 4ea||K|| tx 



QM > 



e 2 1 + 4ea||K| 



Jn e l-4ea|l«lloo7n 



\tp\ 2 dsdt. 
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If we apply this inequality on ip 1 - and if we realize that e 2 J Q |V /, 0~ L | 2 dsdt > e 2 E2\\ip ± \\ 2 , where E2 
is the second eigenvalue of the transverse Laplacian — A^,, we get 



lAr./.i. 



LV J - e 2 V 1 + 4ea||«|| 00 ^ H ' 1 -4ea||«:| 



■* 1 "" + ?*ITe|&" v ^'"- 



Here /3 is a real parameter and from the relation E2 > -Ei it follows that (3 can be chosen in such way that 
for small enough e both coefficients in front of ||'0~ L || 2 an d II^V" 1 "!! 2 are positive. We define a constant 
Cj_ such that the minimum of these two coefficient is equal to ( ^r) which yields 

C± 



ll^ll < e^ v /g^j < sC ± ^QM, 

Using similar ideas, we would get also 



IVVl < e^ X /Q £ A [^] < Cx^jQM- 



ll^ll < sC ± jQ%[4>], 



<eC ± JQ, 



(5.25) 



(5.26) 



(for simplicity we put here the same constant C± even though the estimate could be somewhat finer). 
Now we can finish the estimate on \q(<f>, tp)\ as 



k(M)\ < (C^ k {5(e)) + C 5 e)^Q^}Q^} 



(5.27) 



where C4 and C5 are constants depending on ||k||(x,, C± and a. 

The Hilbert space decomposition will be needed also in the case of last two estimates which are 



technically most difficult. The difference of terms on the second lines of (5.19) and (5.20) reads 

h 



m((j),il)) 



1 



1 V'f V> ds dt - 



Ei 



— 1 (f>ip ds dt. 



Applying formula (5.5) on <f> and ipj we can divide m((^,ip) into four terms. The term m(4>iJ'i,if>iJ'i) 
is integrated by parts with respect to the transverse variable t twice; as usual, we derive the following 
formula for <fi,ip £ Co°(f2o) and we can extend it to W ' (CIq) by density: 

1 f ( h \ r , „ r .,„ „ „, , , 1 f . / h 



m(4>iJi,ipiJi) 



£ ./fir, 



£- - 1 J M1J1 [-A'Ji - E1J1] ds dt - 



^ ./fin 



1 Mi J* dsdt. 



Here the first term vanishes, hence we can estimate 



\m{^iJi^iJi)\< -» 



1 



«=-2 , 
£ -'fio 



-^ 



{fallal J? ds dt 

< 6||«|U / (|fei - fcil + |*a - k" 2 \) \<h\\ik\da + 48e||«||^ /" |0i||Vi|^ 

<12-s/3||K|| O o<rfc(tf(e))||^i||wi.'(/)||^i||i + 48e||K||2 o ||0i||i||^i||j. 
Similarly the terms rn{4>iJi,'4> ) and m(cf>' L , ipiJi) are integrated by parts once to get 

M^iJi,^)! < 8v/3^i (fr fc (<5(£))||0i|| w i,2 (/) + ea || K || oo ||0 1 || / ) IH, 

|™(^iJi,^)| < 8v^ ((r fe (%))|^i|| wl ,=(r)+ea||«|| 00 ||^i||/) JU. 
Finally, the estimate on m(<^ ,^ ) is straightforward, 



(5.28) 
(5.29) 



m(^,V )| < 16ea||K|| 



vv^n || vv 



'„/,-L| 



£1 



l^ll IIV^II 



and (5.26) 



where due to ( 5.25 ) and ( 5.26 1 the term in bracket is bounded. Summing up, we get due to relations ( 5.25 ) 



|m(0,VOI < (Cea k (8(e)) + C 7 e) <J Q^QM , 



(5.30) 

where the constants C§, C7 again depend on ||k||oo, Cj_i a an d in addition on E\. 

Using the Hilbert space decomposition we have shown that in consequence of assumption (2.7) (which 
is automatically satisfied for / bounded), the terms \q(<f>, tp)\ and \m((j>, ip)\ tend to zero in the limit e — > 0. 
In the next section we show why also the assumption (2.8| is needed. 
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5.5.2 Estimates by ag(6(e)) 



Finally, the terms on the first lines of (5.19) and (5.20) are estimated, i.e. we examine the sesquilinear 
form 



KM) 



-^-{d s + od a )^{d s + ed a )^dsdt- f d s 



s ip dsdt — C u 



a2 



<frip ds dt 



(5.31) 



(recall C w = \\d a JiW^i^))- We again have to decompose the functions ip and (j> using (5.5) which leads 



to long expressions where one part of the terms subtracts and other part of terms vanish when e 
due to relations (5.25) and (|5.26|. However, also the following, problematic term occurs 



9d a Ji (i>id s 4> x + <pid 8 il} x ) dsdt. 



(5.32) 



This term is expected to vanish in the limit as e — > due to appearance of (j> ± and ip- 1 . However, 
we do not have any suitable estimate on the longitudinal derivatives of the functions (f>' L and ip , i.e. no 
vanishing control over \\d 3 (f) || and ||<9 S - 0~ L ||. For this reason, we would like to perform an integration by 
parts with respect to the variable s, which requires 9 to be diffcrentiable. To avoid this extra assumption, 
we mollify also 9 using the Stcklov approximation 



9 e (s) 



1 

W) 



s+ 



;i(c) 



5(e) 



9(0 dC 



Here again e t— > 6(e) is a continuous function which vanishes when e — > 0. Then we rewrite the first term 

9 s d a J^id s 4) X dsdt 



in (5.32) in the following way and the integration by parts is justified: 
/ 9d a J 1 i> 1 d s 4>dsdt= [ (d-d^daJ^ds^dsdt^ 



n 

_L 



- 9 £ )d a J 1 ip 1 d s (j)- L dsdt+ / (9 e ) d a Jitpi^ dsdt+ / 9 £ d a Ji(9 s ^i)^ ds dt , (5.33) 



and similarly for the second term in (5.32). In (5.33) the first term will be estimated by <jg(5(s)) and tends 



to zero due to Proposition 4.3 the other terms tend to zero due to (|5.25|, (|5.26) and the boundedness 

(5.34) 



of 9. After tedious computations we get the final formula 

\l(<f>,il>)\ < (c 8 ag(S(e)) + C 9 e) y/Q^QM 
where the constants Cg, and Cg depend on the quantities ||k||ooj C±, a, E\ and also C w , ||^||oo- 

5.6 Summary 

Putting all the estimates (5.23), (5.27), ( |5.30| ) and (5.34) together, we get 

cr(e) := (d + C 2 ) e (\\kt\U + ||fcf |U) + (C 3 + C 5 + C 7 + C 9 ) e+(C 4 + C 6 ) a k (5(e)) + C 8 a d (S(e)) (5.35) 



which tends to zero in the limit e — > due to assumptions of Theorem |2.1| and the proof of this theorem 
is in fact completed. 



In more details, if we implement (5.35) into (5.18) and if we find appropriate constant C as the 

suit 

Let 



of Lemma 



5.2 



maximum of a ll constants involved, we get the statement of Lemma 5.1 In combination with t he r esult 
and relation (5.8) we can set C := C\ [C + C) to get the statement of Theorem 2.1 



us note that the constant C is thus function of A, ||k||oo, C_l, a, -Ed, C u and 



6 Conclusion 



The objective of this paper was to establish the effective Hamiltonian approximation ( 1.1 1 in the norm- 



resolvent sense and under minimal regularity assumptions about the waveguide. Our main result is 



summarized in Theorem 2.1 Let us discuss its assumptions and links with previous results here. 
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6.1 Comparison with previous results 

As mentioned in Section [l.l| the norm resolvent convergence of Theorem |2.1| was proved previously under 
sufficiently regular assumptions about T and 9. Let us show that our conclusions correspond to these 
results. 

In case when k\ , &2 and 9 arc Lipshitz continuous with Lipshitz constants L^ x , L^ 2 and Lg , respec- 
tively, then Hfcfjjoo < I^ku j with jj, = 1, 2, for any choice of 5(e). Consequently, we can abandon the second 
condition in (2.9). Indeed, as we explained in Sectional this assumption was needed to ensure that the 

is not bounded, on the other hand, for Lipshitz continuous fc M 
Hence we can simply cho ose 5(e) — 5(e) = e, 



quantity e||fc„||oo tends to zero if ||fcfj 



this is not the case (cf also the remarks below (5.23)) 

then (Jk(5(e)) oc (L^ + Lk 2 ) e and <jg(5(e)) ex Lae, and similarly as in Theorem 



2.1 



wc find 



U(-A% - e- 2 E 1 - i)- 1 !/- 1 - (ffeff - i)- 1 J 



B(H ) 



<Ce(l + L k 



Lk? + La 



In this way we get the £-type decay rate well known from the other papers (see, e.g., |11]). 

Of course, in case of k\,ki and 9 diffcrcntiable with bounded derivative, Lk t ,Lk 2 and La can be 
replaced by ||fti||<x,, HMoo and l^'Hoo, respectively. 

On the other hand, Theorem |2.1| covers much wider class of curves than previous papers. It is 
reasonable to expect a worse decay rate if the functions fci, &2 and 9 are not differentiable. Then a bound 
on the decay rate can be obtained by optimizing the choice of 5(e), as the following example shows. 



Example 6.1. Let 




s G (2n, In + 1) , 
s € (2n+ l,2ra + 2). 



n G Z . 

Vset. 



The corresponding curve is lying in a plane and is formed by arcs of circle with radius 1 whose center 
is in one half-plane for s G (2ri, In + 1) and in the other half-plane for s G (2n + 1, 2n + 2) (ci the left 
hand side of Figure \M for a part of this infinite curve). Then the best possible estimate reads 



l*fl 



< 



2||M 
6(e) 



2 

W) 



hence on the right hand side of (2.101 the term proportional to eS(e)^ 1 occurs. At the same time, for 
this choice of k\,k,2 it holds o~k(6(e)) ex ^J 5(e). For simplicity we will look for the optimal function 5(e) 
in the class of polynomials, here the most convenient choice is 5(e) — e 2 ' 3 since then the term on the 
right hand side of (2.10) is proportional to e 1 ' 3 (for suitable 9). 



6.2 Optimality of our assumptions 



of the conditions under that our Theorem |2T] holds. 

•2,oo, 



Let us now turn to the optimalitv 

Condition (i) of Assumption 111 requires V G W^°{J; M 3 ), which seems to be the minimal condition to 
guarantee that a (weakly) differentiable moving frame, necessary for the definition of a simultaneously 
twisted and bent tube along the curve, exists. At the same time, the boundedness of curvature k is neces- 



sary to consider the waveguide even as an abstract Riemannian manifold (cf Section 3.2). We therefore 
consider these hypotheses as the natural ones. r\ The same concerns the injectivity assumption (iii) of 
Assumption 111 if we want to interpret the waveguide as a genuine physical device embedded M 3 . But our 
results hold without this last assumption, as pointed out in Remarks |3.3| and |3.4| 

On the other hand, the global boundedness of 9 from Assumption 111 (ii) is not necessary for the 
definition of a non-self-intersecting waveguide. In Figure HI we present an example of an infinite twisted 
waveguide with elliptical cross-section such that 9(s) tends to infinity as \s\ — > oo. It can be introduced 
and handled by the methods of Section|3]without problems. However, the form domain of the transformed 
Laplacian H 6 will not coincide with the Sobolev space W ' (Oq) (i.e. (3.18 1 will not hold) and an extensive 
modification of the present strategy to get the operator limit JO]) would be required. 

2 The case of 'broken-line' waveguide or, more generally, the question of shrinking tubular ncghbourhoods of graphs do 
not fit in the present setting. We refer to recent works of Grieser 1141 1151 for results and references in this field. 
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Figure 4: A waveguide of elliptical cross-section along a straight line whose twisting diverges at infinity. 
The embedded cylindrical channel is responsible for the existence of essential spectrum. 



It is also important to emphasize that the unboundedness of 9 may lead to pathological spectral 
properties. Indeed, the condition 0(s) — > oo as \s\ — > oo implies that H c g has purely discrete spectrum, 
while a css (—A D e ) ^ in general. Actually, the latter happens whenever the cross-section lo contains the 
origin of R 2 (as in Figure H), so that there is an infinite cylindrical channel in fl £ leading to scattering 
waves. It does not contradict the validity of the effective Hamiltonian approximation in principle, since 
the threshold of the essential spectrum of — A^," — e~ 2 Ei tends to infinity in the limit as e — >■ 0, but the 
usefulness of the approximation becomes doubtful. 

We admit that Assumption[2]seems unnatural and it is true that it comes from our technical procedure 
of mollifying the curvature functions fei, &2 an d 0- Although it covers a wide class of waveguides and, in 
particular, all the previously known results, there are still reasonable situations for which Assumption [2] 
does not hold, as the following counterexample shows. 



Example 6.2. Let us define the curve T° 



by giving its curvatures: 






1, S e((n-1 + |K(«-1 + ^) 



•1. •€ (n-l + ^K («-! + ¥) 



In 

2k+2\ 



n€N, fc = 0,l. 



k i r 



2ji 



(n-l)Tt 



Figure 5: The counterexample curve. 

The graph of the function fcj sc is given in Figure [3] as well as T° sc itself. This curve lies in a plane 
and consists of arcs of circle of radius 1 which are shorter and shorter as s grows. For s — > oo ; this curve 
looks like a straight line, however, the curvature is still nonzero. It is possible to show that for all e > 
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there exists jiq G R such that 

sup / \k° sc (s)-k° sc (s + Ti)\ 2 ds = 4ir, VneN, n > n , 

|t 7 |<^-"'('*-i)* 

and i/iis /10/ds true /or any choice of function S. Hence the curvature k° sc does not satisfy Assumption^ 

and our Theorem \2.1\ does not apply. 

To get at least some information about the dynamics in such a pathological quantum waveguide, we 

n osc 
examine the spectrum of the Dirichlet Laplacian — A D e in a untwisted tube with rectangular cross-section 

constructed along the curve T osc . It is possible to show that for any positive e (such that the tube does 

not overlap itself) 

inf a- oss ( -A D ' ^ ) > 0. 



On the other hand, assuming that the effective dynamics is governed by (2.3), we have 

1 



c(il eff ) = er css (H off ) =a[ - A* - - 



-, oo 
4' 



That it, [— \,0] belongs to the essential spectrum of H c g, while the threshold of the essential spectrum of 
the three-dimensional renormalized Hamiltonian is non-negative for every positive e. 

Again, this pathological spectral behavior does not necessarily imply that the norm-resolvent con- 
vergence of — A D C — £~ 2 E\ to H B fi does not hold. As a matter of fact, it is possible to show that 
— Aj-f — £~ 2 E\ possesses an infinite number of negative eigenvalues, hence it may happen that these 
eigenvalues cover the whole interval [— |,0] in the limit as e — > 0. 



In any case, we would like to emphasize that our Theorem |2.1| represents the first norm-resolvent 
convergence result for unbounded waveguides in the full setting of bending and twisting . The question 
of optimality of Assumption [2] in the unbounded case remains open. 

6.3 Two-dimensional waveguides 

The methods of the present paper also enable one to improve the known results [HI E] about the 
effective Hamiltonian approximation in strip-like neighbourhoods of plane curves. The norm-resolvent 



convergence in the two-dimensional case does not follow directly from our three-dimensional Theorem 2.1 
but it can be established exactly in the same way. The proof is in fact much simpler because there is 
just one curvature function, the Frenet frame always exists (it coincides with a relatively parallel frame) 
and there is no twisting (codimension of the reference curve is one). Here we therefore present just the 
ultimate result without proof. The interested reader who is not willing to adapt the present proof to the 
two-dimensional case himself/herself is referred to [21] . 

Let F : / —$■ M 2 be a unit-speed C -smooth curve, where / is an arbitrary open interval (finite, 
semi-infinite, infinite). The vector fields T := T and N := (— T 2 ^ 1 ) form a positively oriented Frenet 
frame of V. We introduce the curvature function k by the Serret-Frenet formula T = —kN. Note that, 
contrary to the three-dimensional case, k is allowed to change sign (and the value of the sign depends 
on the parametrization) 



In analogy with (|3.6[) and (3.7), the two-dimensional waveguide fi £ is introduced as the image (3.6) 



where the mapping C is given now by 

C(s,t) :=r(s)+etN(s) 



with (s,t) E I x (— 1, 1) (see FigureK3j). The unitary transforms ( |2.1[ ) and (2.4) should be replaced by 

U t : L 2 (n e ) -> L 2 {n ,eh{s,t)dsdt) : {V> ^ ?A ° ^l , 
U 2 : L 2 (n ,e 2 h(s,t)dsdt) -> L 2 (fl ) : {ip H- Vehip} , 

where h(s, t) := 1 — en(s)t, and we again define U := U2U1. The latter enables one to approximate in the 
limit as e — > the Dirichlet Laplacian in Q £ by the well known one-dimensional effective Hamiltonian 

k 2 

Heft '■= —A n r ■ 



2G 



N(s) 




t=l 



Figure 6: A two-dimensional quantum waveguide. 



The main idea of the proof again consists in replacing Ui by (4.4) using the mollification k £ defined 



analogously to (2.5) 



The two-dimensional version of Theorem 12. f I reads as follows. 
Theorem 6.3. Let the following assumptions hold true: 
(i) re Wf ' c °°(/;IR 2 ) and k<=L°°(I). 
(ii) Q e does not overlap itself (i.e. C is injective) for small enough e. 



(hi) \im a K (5(e)) = for some positive continuous function 5 satisfying (4.1), where a K is defined 



by f^ef. 

Then there exist positive constants £q and C such that for all e < Eq, 



U{-A% - e- 2 E x - i)- l U- 1 - (H cS - i)- 1 © 0- 1 



< 



B(L 2 (fi )) 



c(e + s\\k E \\+o- K (6(e))) 



where 1 " denotes the zero operator on the orthogonal complement of the span of {(p C3> 3\ \ <-P € L 2 (I)} 
andU = U 2 U 1 . 



6.4 Different boundary conditions 



It is well known that the structure of the effective Hamiltonian in the limit ( 1.1 ) is a consequence of the 



choice of Dirichlet boundary conditions on the 'lateral boundary' / x du). Indeed, there is no geometric 
potential for Neumann boundary conditions 22J and the limit can have a completely different nature 
if one considers a combination of Dirichlet and Neuman boundary conditions |18j . On the other hand, 
our choice of Dirichlet boundary conditions on (dl) x uj is not essential. In the same manner, we could 
impose any other type of boundary conditions (Neumann, Robin, periodic, etc), which would lead to 
an analogue of Theorem |2.1[ with the boundary conditions of H c s changed accordingly. In particular, 
the choice of periodic boundary conditions enable us to cover the case of tubes about closed (compact) 
curves. 



6.5 Non-thin waveguides 

Finally, let us mention that the tricks of the present paper how to deal with quantum waveguides under 
mild regularity assumptions do not restrict to the effective Hamiltonian approximation. For instance, 
the usage of the relatively parallel frame instead of the Frenet frame enables one to extend some spectral 
results for — A^ (with e not necessarily small) to the general case of waveguides along merely twice 
differentiable curves with possibly vanishing curvature. In particular, we have in mind the classical 
results about the curvature- induced bound states 9, 5J and the recent ones about Hardy- type inequalities 
due to twisting [TU1 IPTI IT!?] . 



27 



Acknowledgement 

We are grateful to Denis Borisov for valuable discussions, in particular for his idea how to significantly 
relax our initial hypotheses about the waveguide regularity This work has been partially supported by 
the Czech Ministry of Education, Youth, and Sports within the project LC06002 and by the GACR 
grant No. P203/11/0701. 

References 



[i] 

[2] 
[3] 

[4] 

[5] 
[6] 

[7] 
[8] 

[9 
[10 

[11 

[12 

[13 
[14 

[15 

[16 



N. I. Akhiezer, Theory of approximation, Frederick Ungar Publishing Co., New York, 1956, 2nd 
augmented Russian edition is available, Moscow, Nauka, 1965. 

R. L. Bishop, There is more than one way to frame a curve, The American Mathematical Monthly 
82 (1975), 246-251. 

D. Borisov and G. Cardone, Complete asymptotic expansions for the eigenvalues of the Dirichlet 
Laplacian in thin three-dimensional rods, ESAIM: Control, Optimisation, and Calculus of Variation, 
to appear. 

G. Bouchitte, M. L. Mascarenhas, and L. Trabucho, On the curvature and torsion effects in one 
dimensional waveguides, ESAIM: Control, Optimisation and Calculus of Variations 13 (2007), 793- 
808. 

B. Chenaud, P. Duclos, P. Freitas, and D. Krejciffk, Geometrically induced discrete spectrum in 
curved tubes, Differential Geom. Appl. 23 (2005), no. 2, 95-105. 

C. R. de Oliveira, Quantum singular operator limits of thin Dirichlet tubes via T -convergence, Rep. 
Math. Phys. 66 (2010), 375-406. 

C. R. de Oliveira and A. A. Verri, On the spectrum and weakly effective operator for Dirichlet 
Laplacian in thin deformed tubes, J. Math. Anal. Appl. 381 (2011), 454-468. 

C. R. de Oliveira and A. A. Verri, On norm resolvent and guadratic form convergences in asymptotic 
thin spatial waveguides, to appear in Proceedings of Spectral Days 2010, Santiago, Chile, in the book 
series "Operator Theory: Advances and Application" , Birkhiiser, Basel. 

P. Duclos and P. Exner, Curvature-induced bound states in guantum waveguides in two and three 
dimensions, Rev. Math. Phys. 7 (1995), 73-102. 

T. Ekholm, H. Kovafik, and D. Krejciffk, A Hardy inequality in twisted waveguides, Arch. Ration. 
Mech. Anal. 188 (2008), 245-264. 

P. Freitas and D. Krejciffk, Location of the nodal set for thin curved tubes, Indiana Univ. Math. J. 
57 (2008), no. 1, 343-376. 

L. Friedlander and M. Solomyak, On the spectrum of the Dirichlet Laplacian in a narrow infinite 
strip, Amcr. Math. Soc. Transl. 225 (2008), 103-116. 

, On the spectrum of the Dirichlet Laplacian in a narrow strip, Israeli Math. J. 170 (2009), 

no. 1, 337-354. 

D. Grieser, Thin tubes in mathematical physics, global analysis and spectral geometry, Analysis on 
Graphs and its Applications, Cambridge, 2007 (P. Exner et al., ed.), Proc. Sympos. Pure Math., 
vol. 77, Amer. Math. Soc, Providence, RI, 2008, pp. 565-593. 

D. Grieser, Spectra of graph neighborhoods and scattering, Proc. London Math. Soc. 97 (2008), 
718-752. 

T. Kato, Perturbation theory for linear operators, Springer- Verlag, Berlin, 1966. 



28 



[17] D. Krejcifik, Twisting versus bending in quantum waveguides, Analysis on Graphs and its Applica- 
tions, Cambridge, 2007 (P. Exner et al., ed.), Proc. Sympos. Pure Math., vol. 77, Amer. Math. Soc., 
Providence, RI, 2008, pp. 617-636. See arXiv:0712.3371v2 [math-ph] (2009) for a corrected version. 

[18] D. Krejcifik, Spectrum of the Laplacian in a narrow curved strip with combined Dirichlet and Neu- 
mann boundary conditions, ESAIM: Control, Optimisation and Calculus of Variations 15 (2009), 
555-568. 

[19] D. Krejcifik and E. Zuazua, The Hardy inequality and the heat equation in twisted tubes, J. Math. 
Pures Appl. 94 (2010), 277-303. 

[20] J. Lampart, S. Teufel, and J. Wachsmuth, Effective Hamiltonians for thin Dirichlet tubes with 
varying cross-section, Mathematical Results in Quantum Physics, Hradec Kralove, 2010, xi+274 p., 
World Scientific, Singapore, 2011, pp. 183-189. 

[21] H. Sedivakova, Quantum waveguides under mild regularity assumptions, Diploma Thesis, Czech 
Technical University, June 2011, Supervisor: D. Krejcifik. The electronic version is available at 
http: //ssmf . f jf i . cvut . cz/studthes/2008/Sedivakova_thesis . pdf . 

[22] M. Schatzman, On the eigenvalues of the Laplace operator on a thin set with Neumann boundary 
conditions, Applicable Anal. 61 (1996), 293-306. 

[23] J. Wachsmuth and S. Teufel, Effective Hamiltonians for constrained quantum systems, preprint on 
arXiv:0907.0351v3 [math-ph] (2009). 

[24] , Constrained quantum systems as an adiabatic problem, Phys. Rev. A 82 (2010), 022112. 



29 



